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Abstract 



Under general conditions, the equation g{x, y) = implicitly defines y locally as 

a function of x. In this short note we study the combinatorial structure underlying 

^ a recently discovered formula for the divided differences of y expressed in terms of 

^ bivariate divided differences of 5, by analyzing the number of terms a„ in this formula. 

0^ The main result describes six equivalent characterizations of the sequence {a„}. 

O 

1 Introduction 

o 

CN The act of taking differences and dividing is a recursive process with a wide range of ap- 

plications in mathematics. In approximation theory, such divided differences can be viewed 
as a discrete analogue of derivatives; see the survey by de Boor [2]. Generalizations of uni- 

^ variate divided differences play a role in the theory of orthogonal polynomials, Macdonald 

^ polynomials in particular, and are used to define recurrences between Schubert polynomials 

in the algebraic geometry of flag varieties [11]. 

Although a classical topic dating back to Newton, interest in divided differences seems 
to have been revived recently. One point of view deflnes multivariate divided differences as 
the leading term of an interpolating polynomial and attempts to carry over properties from 
the univariate case [15, 17]. In addition, several formulas from (differential) calculus have 
been generalized to a 'divided difference calculus', including several new chain rules [6, 7, 21] 
akin to Faa di Bruno's formula, and expressions for divided differences of inverse functions 
[8] and implicit functions [12, 14]. 
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It is the goal of this paper to study the combinatorial structure underlying a formula by 
Muntingh and Floater [12], restated below as (5), for the divided differences of a function 
y implicitly defined by a relation g. More precisely, for some open intervals [/, C M, let 
y : U — y be a function that is implicitly defined by a function g : U x V — > M via 

g{x,y{x))=0, ^{x,y{x))^0 y x e U. (1) 

Using a chain rule for divided differences, this relation induces relations between the divided 
differences of y and those of g and in particular yields (5). In this paper we study the number 
of terms in this formula. 

In the following section, we recall the definition of a divided difference, a recurrence 
relation for the divided differences of y, and the explicit formula (5). Next, in Section 3, we 
first derive a generating function from the recurrence relation and then prove six equivalent 
characterizations of the sequence {a„}. We end with an analysis of the asymptotic behaviour 
of {an}. 



2 Divided Differences 

Let [xo, . . . , Xn]f denote the divided difference [of order n] of a function / : (a, b) — > M at 
the distinct points xq, . . . , a;„ G (a, b), which is recursively defined by [xo]f := /(xq) and 

[Xo, . . . , X„]/ = [''^^■■■^^n]f-[Xo,..-.Xn^,]f -^^^q 

Xji Xq 

For given indices io,ii, . . . ,ik satisfying io < ii < ■ ■ ■ < it, v/e shall shorten notation to 
[ioii---ik]f ■■= [xio,Xi,,...XiJ/. 

The above definitions generalize to bivariate divided differences as follows. This time, let 
/ : U — y M be defined on some rectangle 

U = (ai,6i) X (a2,&2) C Ml 

Suppose we are given integers m,n > 0, distinct points xq, . . . ,Xm G (ai,6i), and distinct 
points yo, . . . ,ym & (02, &2)- The Cartesian product 

{xo, . . . ,Xm} X {yo, ...,yn} 

defines a rectangular grid of points in U. The [bivariate] divided difference of / at this grid, 
denoted by 

[xQ,...,Xm]yo,---,yn]f, (2) 

can be defined recursively as follows. U m = n = 0, the grid consists of only one point 
{xo,yo), and we define [xo;i/o]/ •= f{xo,yo) as the value of / at this point. In case m > 0, 
we can define (2) as 

[Xi, . . . , Xm] • • • 5 Vnlf — [Xq, ■ ■ ■ , Xm-l', UOi ■ ■ ■ i Vnlf 
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Figure 1: The dissections of a convex polygon with vertices labeled 0, 1, 2, 3,4. The bottom- 
right dissection has faces (0,1,4), (1,2,4), and (2,3,4). 



or if n > 0, as 



[xo, . . . , yi, . . . , yn]f - [xo, . . . , x^; ?/o, • • • , l/n-i]/ 



Vn - 2/0 

If both m > and n > the divided difference (2) is uniquely defined by either recurrence 
relation. Similarly to the univariate case, we shorten the notation for bivariate divided 
differences to 

[ioii ■ ■ - is] joji ■ ■ ■ jt]f '■= [xio,Xi-i^, . . . , Xjj yjQ,yj-^, . . . , 

In what follows, assume that y and g are related by (1). Applying a chain rule for divided 
differences, one can derive [12] 

' ' |0; Oils' 
and for any n > 2 the recurrence relation 

[oi...„b=-5: E ^^T^tF^niw.-.+i)---*. (4) 

k=2 0=io<--<ik=n s=0 ^ ' l=s+l 

s=is—io 

The third summation might look a bit mysterious. It is a concise way of describing all s for 
which the increasing sequence = io < "^i < ■ ■ ■ < "^fc = starts with at least s steps of 1, 
i.e., h - io = ■ ■ ■ = is - is-i = 1- 

For n > 2, this recurrence relation was "solved" to yield the explicit formula (5). In order 
to describe this formula, we introduce some notation for dissections of a convex polygon. 
With a sequence of labels 0, 1, . . . , n we associate the ordered vertices of a convex polygon. 
A dissection of a convex polygon is the result of connecting any pairs of nonadjacent vertices 
with straight line segments, none of which intersect. We denote the set of all such dissections 
of a polygon with vertices 0,1, ... ,n by V{0, 1, . . . ,n). The vertices 0,1, ... ,n and line 
segments (between either adjacent or nonadjacent vertices) form the vertices and edges of 
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n 


: 


1 2 


3 


4 


5 


6 7 


8 


9 


10 


11 


On 


: 1 


1 3 


13 


71 


441 


2955 20805 


151695 


1135345 


8671763 


67320573 



Table 1: The number of terms a„ in the right hand side of Equation (5) for small values of 
n, with the additional convention 



a plane graph. As such, every dissection tt G 7^(0, 1, . . . , n) is described by its set -F(vr) 
of [oriented] faces, which does not include the unbounded face. Each face / G -F(vr) is 
represented by a subsequence / = (io, ii, ■ ■ ■ ,ik) of the sequence (0, 1, ... , n) of length at 
least three. 

See Figure 1 for the 11 dissections of a convex polygon with vertices 0,1,2,3,4. For 
general n, this number is given by the little Schroder number [19] (or super-Catalan number) 

1 ^ /n - 1\ /ri - 1 + 

V k )\ k-1 

k=i ^ ^ ^ 

which form sequence A001003 in The On-Line Encyclopedia for Integer Sequences (OEIS) [18] . 
The main result of [12] is the formula 

[01- ■■n]y= y T] y -i^^^i^^i^v^ n (.ihiiiiMay 

S=is—io i;— i;_l=l 

which expresses divided differences of y solely in terms of bivariate divided differences of g. 
In a follow-up paper, this formula was generalized to a formula for divided differences of 
multivariate implicit functions [13, 14]. 



3 The Number of Terms 

For any implicit divided differences [01 ■ ■ - n]?/ of any order n > 2, we wish to analyze the 
number of terms a„ in the right hand side of Equation (5). Equation (3) yields the natural 
extension ai = 1. With the additional convention oq = 1, we can collect these integers 
into a sequence {a„}. The sequences in this paper will always start at n = 0, unless stated 
otherwise. Alternatively, we can encode the integers by an [ordinary] generating function, 
which is the formal power series Yl'^=o ^nX^- Table 1 lists for small values of n, and more 
terms can be found as sequence A162326 in The OEIS. 

To analyze {a^}, we rewrite (4) in terms of the more standard notion of a composition. 
Let N = {1,2, . . .} denote the set of positive integers. A composition [of n of length k] is 
any sequence (ji, . . . ,jk) G satisfying ji + ■ ■ ■ + jk = n. The positive integers ji, . . . ,jk 
are called the parts of the composition. Clearly these correspond to increasing sequences 
= io < ■ ■ ■ < ik = n via ji = ii — ii-i for I = 1, . . . , k. Thus (4) is equivalent to 

n k 

[01...n]y = -J2 E E (6) 

k=2 compositions s=0 

ji+-+jk=n s=ji+-+js 
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[01- ■■s;{s+js+i)-- -nig -|-|- . , n 
[(^1 + ■ ■ ■ + ^'-) ■ ■ ■ Oi + ■ ■ ■ + ^0] 2/. 

Proposition 1. The sequence {an} has associated generating function 



oo 



n=0 



Proof. Replacing, in Equation (6), the divided differences of y of order n by tlie number of 
terms an yields 



k 



k=2 compositions s=0 1=1 
h+-+3k=n s=h+-+js 

The third sum makes sure that any composition starting with precisely s ones will be counted 
s + 1 times. We can, therefore, replace this sum by a sum over the number t of initial parts 
in the composition that are one. Using that ai = 1, we find 

n k—1 k—t 

an = l + J2 Yl Y n ^i'' ^ - 2' 

k=2 t=0 compositions 1=1 

ji-i ^ik-t=n-t 

where the first term, 1, is included to make sure that the composition n = l + l + -- - + lis 
counted n + 1 times. The generating function G will therefore satisfy 



k=2 t=0 
^ oo k 

Y^ + iG-^)YY^'(^G-ir^ 



1 — X 



k=l t=0 

1 1 



1-X2-G 



- 1 



Multiplying by (1 — x)(2 — G) and rewriting gives 

2G^ -5G + 2 + ^— = 0. (8) 
1 — X 

Since G{0) = 1, we find that G is as in the Proposition. □ 

Up to a constant and a rescaling, the generating function G{x) is exactly the magnetiza- 
tion Moo(s) for the XYZ spin chain along a special line of couplings [4]. 

Before we state the next proposition, we need to recall some standard notions. The 

Gaussian hypergeometric series [with parameters a, b, c] is defined by 

IT. / r A {a)n{b)n 2:" 

2Fiia,b;c;z) = }^—-—-, 

n=0 ^^^"^ 



where the Pochhammer symbol (■)„ is defined by 

1, ifn = 0; 

a(a + 1) ■ ■ ■ (a + n — 1), if > 0. 

Hypergeometric series with adjacent parameter values are related [1, Section 15.2.11] by 
Gauss' contiguous relation 

(c - b)2Fi{a, b-l;c]z) + {2b-c+{a- b)z)2Fi{a, b; c; z) + b{z - 1)2^1(0, 6+ 1; c; z) = 0. (9) 

Moreover, as explained by Roy [16], sums of products of binomial coefficients can often be 
reduced to hypergeometric series, with the help of elementary identities like 



Given a sequence {a„}, define its binomial transform by 

6„ = ^Qafc, n>0. 



It is well known that if A{x) is the generating function of an integer sequence {a„}, then 

1 X 



1 — X \1 — a; 

is the generating function of its binomial transform. Although sometimes proved in the case 
that A{x) is analytic [3, Appendix], the statement holds for any formal power series A{x) [9]. 
For any n > 0, let 

„ 1 /2n^ 



n + 1 \ n 

denote the n-th Catalan number, which form sequence A000108 in The OEIS. It is well known 
[10] that the sequence {C„} has generating function 

, ^ 1 — \/l — Ax , , 

Fix) := —^^ . (11) 

By means of the binomial transform, we shall relate {a„} to {2"C„}, which form sequence 
A151374 in The OEIS. 

Proposition 2. With the convention 1, the following statements are equivalent 

and hold. 

(a) The number of terms in the r.h.s. of (5) is an for any n > 2; 

(b) The sequence {an} has generating function given by (7); 
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(c) The integers an satisfy the quadratic recurrence relation 



n-l 



Ifi 1 -|- 2 ^ ^ aman—mi ^ 2, 



(d) The integers an satisfy the linear recurrence relation 

nan = (-14 + 10n)a„_i + (18 - 9r;,)a„_2, n>2] 

(e) One has an = 2-^1 ( 2' ~ ^' ~l ^ ~ 



n>0' 



(f) The sequence {a„+i}„>o is the binomial transform o/|2"'C„} 

Proof. Clearly (a) holds by definition of the {an}, and (a) <^=^ (b) by Proposition 1. 

"(b) <^=^ (c)": Given that oq = 1 and Oi = 1, the sequence {an} has generating function 
G given by (7) if and only if (8) holds, i.e., if and only if 

(00 \ ^ 00 00 

5]]a„x" -5^a„a;" + 2 + ^a;" 
ra=0 / n=0 n=0 

00 / n \ 

= ^ 1 - SOn + 2 ^ aman-m ^ + 2 
n=0 \ m=0 
oo / n—l 

= ^ l-a„ + 2^ anittn- 



n=2 \ m=l 



which happens precisely when {a„} satisfies (c). 

"(b) =^ (d)": In the differential ring of formal power series with derivation 

the generating function G satisfies the differential equation 

= (1 - lOx + 9x^)G' + 4^-5, 

which can be written in terms of the {an} as 

= (—5 + 4ao + fli) ■ 1 + (2a2 — 6ai) ■ x 

oo 

+ (^(n + l)a„+i + (4 - 10n)a„ + 9(n - l)a, 

n=2 



1 ■ X 



Comparing coefficients yields (d). 

"(d) ^> (e)": Substituting a = |, 6 = 2 — n, c = 2, 2 = — 8 in Equation (9), one finds 
that 2F1 (|, 1 — n;2; —8) satisfies the same recurrence relation as a^, and must therefore be 
equal to a„. 

"(e) =^ (f)": The binomial transform of {2"C„} is given by 
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Using the identities (10) with a = 1, one finds 

^ (2A;)! n! 1 



k\ {n-k)\ {k + l)\ k\ 

_ ^ CM-^h (-8)^ 

^ (2)fc A;! 
= Q,-n;2;-8 

for any n > 0. 

"(f) =^ (b)": Since {C„} has generating function F(x) given by (11), the sequence 
{2"C„} has generating function 



and its binomial transform {an+i}n>o has generating function 



8x 1 /l-9x 



n=0 

We conclude that {a„} has generating function 



Ax Ax 



oo 



1 / l-9x \ I 

1 - V 13^ \ 5 1 /I - 9x 



n=0 



We end this paper with a remark on the asymptotic behaviour of the sequence {a„}. 
Flajolet and Sedgewick [5] describe how the asymptotic properties of a sequence can be 
deduced from the singularities of its generating function, i.e., points where the function ceases 
to be analytic. The location p of the singularity closest to the origin dictates an exponential 
growth p"", while the singularity type determines a subexponential growth factor C{n). 
Using the dictionary in Figure VI. 5 in [5], the expansion 

around the singularity at p = | yields the asymptotic formula 

C{n) ■ 9", C{n) = -^^n-^/^ + 0{n~'/^), (12) 

16v27r 
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Figure 2: The relative error 1 — C{n) ■ 9"/a„ plotted against n. 



for the coefficients a„. Using Proposition 2(d) to quickly calculate the ffist 1000 entries in 
{an}, we compare (12) to by plotting the relative error 

^ _ C{n) ■ 9" 

in Figure 2. 
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